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ABSTRACT 

The nearby M-dwarf AP Col was recently identified by Ricdcl ct al. 2011 as a pre-main-sequence 
star (age 12 - 50 Myr) situated only 8.4 pc from the Sun. The combination of its youth, distance, and 
intrinsically low luminosity make it an ideal target to search for extrasolar planets using direct imaging. 
We report deep adaptive optics observations of AP Col taken with VLT/NACO and Keck/NIRC2 in 
the L-band. Using aggressive speckle suppression and background subtraction techniques, we are 
able to rule out companions with mass m > 0.5 - lMj up for projected separations a > 4.5 AU, and 
m > 2 Mj up for projected separations as small as 3 AU, assuming an age of 40 Myr using the COND 
theoretical evolutionary models. Using a different set of models the mass limits increase by a factor of 
>2. The observations presented here are the deepest mass-sensitivity limits yet achieved within 20 AU 
on a star with direct imaging. While Doppler radial velocity surveys have shown that Jovian bodies 
with close-in orbits are rare around M-dwarfs, gravitational microlensing studies predict that 17~tg% of 
these stars host massive planets with orbital separations of 1-10 AU. Sensitive high-contrast imaging 
observations, like those presented here, will help to validate results from complementary detection 
techniques by determining the frequency of gas giant planets on wide orbits around M-dwarfs. 

Subject headings: stars: formation — planets and satellites: formation — planets and satellites: 
detection — stars: individual (AP Col) 



1. INTRODUCTION 

The high-contrast imaging technique has provided 
some remarkable extrasolar planet discoveries over the 
past few years (e.g., Chauvin et al. 2005; Marois et al. 
2008, 2010; Lagrange et al. 2010). Direct detection allows 
us not only to infer orbital separations and eventually 
determine orbits (Crepp et al. 2011), but also to study 
extrasolar planets in detail through multi-color photom- 
etry and spectroscopy (Mohanty et al. 2007; Quanz et al. 
2010; Hinz et al. 2010; Janson et al. 2010; Bowler et al. 
2010). Investigation of the small sample of directly im- 
aged planets known to date has already revealed unex- 
pected results, including the suggested presence of thick 
cloud layers (Skemer et al. 2011), enhanced metallicity, 
and indications of non-equilibrium chemistry (Barman 
et al. 2011). Such observations present a challenge to 
theoretical atmospheric models (Skemer et al. 2012). 

Dedicated imaging surveys have now searched host 
stars of all spectral-types, with an initial emphasis on 
FGK stars (Masciadri et al. 2005; Biller et al. 2007; 
Kasper et al. 2007; Lafreniere et al. 2007a; Chauvin et al. 
2010; Heinze et al. 2010). The fact that the first hand- 
ful of directly imaged planets have been found orbiting 
A-stars provides an unambiguous clue that the efficiency 
of the planet formation process depends upon host star 
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mass. Indeed, extrapolation of the Doppler radial veloc- 
ity planet population to separations accessible to high- 
contrast imaging instruments predicts that A-type stars 
are ideal targets, even though they are intrinsically bright 
(Crepp & Johnson 2011). Massive stars have a high gas 
giant planet occurrence rate, and a propensity to form 
more massive planets with long orbital periods (John- 
son et al. 2007). While the same line of reasoning sug- 
gests a relative paucity of giant planets orbiting M-stars 
(Kennedy & Kenyon 2008), initial results from gravita- 
tional microlensing surveys suggest a significantly higher 
frequency of planets located beyond the snow- line (Gould 
et al. 2010; Cassan et al. 2012) compared to rates of close- 
in planets inferred from Doppler measurements (Johnson 
et al. 2010; Bonfils et al. 2011). 

High-contrast imaging observations will soon be able 
to rectify any discrepancies found by indirect planet de- 
tection techniques. Planned instruments and new adap- 
tive optics (AO) systems at the VLT (Beuzit et al. 2006), 
Gemini South (Macintosh et al. 2007), Palomar (Hinkley 
et al. 2011), the LBT (Esposito et al. 2012), and Subaru 
(Tamura 2009) will generate unprecedented contrast lev- 
els, providing direct access to young planets with mass 
m ~ 0.5Afj up . It is imperative to reach this level of sensi- 
tivity because the "bottom-to-top" formation paradigm 
of core-accretion (Pollack et al. 1996) predicts an overall 
planet occurrence rate that rises steadily with decreasing 
planet mass. In other words, access to sub- Jovian-mass 
bodies significantly enhances the prospects for detecting 
an extrasolar planet around any type of star. 

Youth and distance are important parameters to con- 
sider when prioritizing high-contrast imaging targets. 
With an accurately measured distance of only 8.39±0.07 
pc and estimated age between 12 - 50 Myr, the nearby 
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M4.5 flare star AP Col (RA: 06 /l 04 m 52M6; DEC: 
-34°33'36".0 (J2000, E2000)) has recently been iden- 
tified as the closest pre-main sequence star to the Sun 
(Riedel et al. 2011). Based on its apparent age and space 
motion, AP Col is most likely a member of the —40 Myr 
old Argus/IC 2391 Association (Riedel et al. 2011). Stel- 
lar companions have been ruled out by lucky imaging ob- 
servations prior to establishing its proximity and youth 
(Bcrgfors et al. 2010). Given its (extremely) convenient 
properties, AP Col is an ideal target to search for gas 
giant planets by means of direct imaging. 

In this paper, we present the results from two different 
sets of deep high-contrast imaging observations taken at 
the VLT and Keck. The observations were carried out 
in the L-band which, compared to shorter near-infrared 
(NIR) wavelengths, provides a higher Strehl ratio and 
more favorable star-planet flux ratio, while also main- 
taining manageable thermal background levels from the 
sky and instrument optics. Despite generating the deep- 
est direct imaging mass-sensitivity yet achieved for or- 
bital separations on the scale of the solar-system (a > 3 
AU), we find no evidence for the existence of gas giant 
planets. We discuss the implications of these observa- 
tions to the extent that they support the notion that 
Jovian bodies appear to be rare around M-dwarfs. 

2. OBSERVATIONS AND DATA REDUCTION 

2.1. VLT/NACO 

Imaging data was acquired on December 4, 2011 UT, 
with the AO-fed, high spatial resolution camera NACO 
mounted on ESO's VLT/UT4 8.2-m telescope at Paranal 
(Lenzen et al. 2003; Rousset et al. 2003). All images were 
unocculted and taken with the L27 camera (plate scale 
~ 27.15 mas pixel -1 ) using the L' filter (A c = 3.8/zm, 
AA = 0.62 /im). We used pupil stabilization mode to 
enable angular differential imaging (ADI) (Marois et al. 
2006). 

To enhance the signal-to-noise ratio (S/N) of potential 
companions, we chose to moderately saturate the stel- 
lar PSF core. Unsaturated images were also acquired to 
calibrate the photometry. Detector reads were recorded 
individually ("cube mode"). The integration time was 
set to 0.2s per read (unsaturated reads had an expo- 
sure time of 0.1s). We obtained 134 raw data cubes for 
AP Col in total. Each cube consisted of 64 reads. Be- 
tween cubes we moved the telescope by «9" on the sky 
executing a 5-point dither pattern to facilitate removal 
of background noise from the sky and instrument op- 
tics. Individual detector reads were checked for open 
loops and poor AO correction. We discarded one com- 
plete cube as well as several individual reads from other 
cubes. The remaining reads were median combined in 
stacks of 32, yielding 265 science images with an effective 
integration time of 32x0.2 s each. We measure the PSF 
FWHM in unsaturated images to be -4.3 px (« 0.12"), 
which is comparable to the theoretical diffraction limit, 
O ss \ C /D m 0.10". 

Our basic data reduction steps (bad pixel correction, 
sky subtraction) are described in Quanz et al. (2010, 
2011). We use the LOCI algorithm to subtract the 
stellar PSF (Lafreniere et al. 2007b). Following the 
same naming convention as in the original paper, we set 
pertinent reduction parameters to the following values: 
FWHM=4.5 px, A,5=0.75, dr=5, AT A =300. 



2.2. Keck/NIRC2 data 

We also observed AP Col from Keck Observatory on 
January 7, 2012 UT. Using the 10-m Keck II telescope 
and AO system (Wizinowich et al. 2004) , we acquired im- 
ages with the Near InfraRed Camera (NIRC2; PI: Keith 
Matthews) in the L' filter (A c = 3.8^m, AA = 0.7 pm). 
A total of 67 images were obtained, each consisting of 
300 reads with 0.106s integration time per read. All im- 
ages were in the linear detector regime and recorded with 
the narrow camera setting which provides a plate scale of 
9.963 mas pixel" 1 (Chez et al. 2008). Like the VLT data, 
we did not use a coronagraph. Vertical angle mode was 
used to allow for ADI operation. The AO system was 
running at a refresh rate of 438 Hz. The FWHM of the 
PSF was 9.1 px (w 0.09"), i.e., close to the theoretical 
diffraction limit of 9 w 0.08". 

Raw images were processed using standard techniques 
to remove dark current, flat-field, replace hot-pixel val- 
ues, and align the target star. In the same manner as 
described in the previous section with the VLT data, we 
used the LOCI algorithm on the 67 images to subtract 
the stellar PSF. The only difference was the value for the 
FWHM which was set to 9 px for the Keck data. Table 1 
shows a summary of relevant observing parameters for 
both the VLT and the Keck datasets. 

3. RESULTS 

We find no evidence for point sources in either the 
VLT or Keck data sets, in raw or processed images. Our 
checks for real companions did initially show potential 
candidates in the VLT data. However, they subsequently 
disappeared when combining only a fraction of the final 
images, and other bright speckles showed up instead. 

3.1. Derivation of contrast curve 

A comparison of the 5-a contrast levels obtained with 
NACO and NIRC2 is shown Figure 1. Contrast is defined 
as: 

AL'(r) = 2.5 • log 10 ^T= W 

5cr(r)/yV2 p 

where is the mean stellar flux per pixel in an aper- 
ture with radius r ap = 3 pixels (VLT/NACO) or r ap = 5 
pixels (Keck/NIRC2), and a is the standard deviation 
of the pixels within centro-symmetric annuli centered 
on the star, which are twice as wide as the aperture 
radius. The VLT/NACO observations provide slightly 
deeper contrast compared to Keck as a result of lower 
background levels, more field rotation, and lower airmass. 
We focus on these data for the remainder of the analysis. 

To derive our final detection limits, we insert fake plan- 
ets of known brightness into individual NACO frames. 
Two planets of the same brightness are inserted at the 
same radial distance (0.3", 0.4", 0.5", 0.6", 0.8", 1.0", 
1.5", and 2.0"), but at opposing sides from the central 
star. The signal-to-noise of the planets in the final LOCI 
processed image is computed via: 

( — \ -^Planet 
^Planet „ / J^T 

where i^pianot is the measured planet flux in an aperture 
with radius r ap = 3 pixels and a is the standard deviation 
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of the pixels within an 7 pixel-wide arc between the two 
fake planets and centered at the same radial distance. 
By inserting two planets simultaneously and having two 
arcs for the noise estimate between the planets, we ob- 
tain 4 different values for (^) Planet - To be conservative, 
we chose the combination yielding the lowest value as our 
reference. For a given radial separation we lowered the 
brightness of the planets in steps of 0.5 magnitudes until 
(jf) Planet < ^" ^ e s ig na l-to- noise of the next brightest 
planet was then used for calibration and the values ex- 
trapolated to a level that corresponds to (^) Planot = 5- 
These values are shown in Figure 2. 

We emphasize that injecting fake planets of various 
brightness and retrieving them is the only way to ob- 
tain an accurate estimate for the achieved detection 
limit. LOCI removes flux from potential companions 
(Lafreniere et al. 2007b) as a function of separation from 
the star, LOCI input parameters, and brightness of the 
companions themselves (Pueyo et al. 2012). There are 
no analytical methods for predicting and correcting for 
flux losses. 

The limiting magnitude in the VLT/NACO data ap- 
proaches L[ imit w 16.7 mag at a separation of « 0.7" , be- 
yond which the sensitivity becomes flat indicating that 
we have reached the floor set by thermal background 
noise (Fig. 2). Sensitivity is dominated by subtraction 
residuals closer to the star where the detection limits are 
governed more-so by residual scattered light. 

3.2. Detection limits using COND models 

To convert our on-sky contrast into mass-sensitivity, 
we must first estimate the V magnitude of AP Col. Us- 
ing color transformations provided by the 2MASS team 5 , 
we estimate AP Col's K magnitude to be K = 6.88±0.05 
mag, based on the apparent K 8 value of K s = 6.87±0.02 
mag (Cutri et al. 2003). A typical M4.5 star has an ap- 
proximate K — L' color of 0.4 mag (Cox 2000). Thus, 
AP Col has an V magnitude of V « 6.48 ± 0.05 mag. 

Figure 2 (left-hand panel) compares our detection lim- 
its with the predicted L' magnitudes for 40 Myr objects 
having masses of 0.5, 1 and 2 Mj up based on the COND 
models (Baraffe et al. 2003). We select this specific age 
because Riedel et al. (2011) derive a most likely age be- 
tween 12 - 50 Myr for AP Col, and membership with the 
~40 Myr old Argus/IC 2391 association seems plausible 
from AP Col's space motion. Our data are sufficiently 
sensitive to detect objects with masses between 0.5 - 1 
Mjup (or greater) for separations >4.5 AU. Inward of 4.5 
AU, as close as ~3 AU, objects with masses between 1-2 
Mjup (or greater) are detectable. Based on these model 
predictions, this is the first time that a young Jupiter 
analog could have been detected by means of direct imag- 
ing around another star. This is also demonstrated in 
Figure 3, where we have inserted two fake planets in the 
raw images at a projected separation of 5.2 AU and re- 
run the data reduction. The planets' brightness corre- 
sponds to a 1 Mj up object with an age of 40 Myr based 
on the COND models. It clearly shows that we would 
have detected such an object in our data. 

To assess the likelihood of imaging a planet during a 
single epoch observation, the following relation can be 

5 http://www.astro.caltech.edu/ jmc/2mass/v3/transformations/ 



used 

P, detect (a) = cos(arcsin(IWA/a)) , (3) 

where Pdotect is the detection probability, a the semi- 
major axis of the planet in AU, and IWA is the inner 
working angle (in AU) where the contrast curves allows 
for the detection of this planet. This relation is only 
valid for circular orbits and assumes uniform sensitivity 
in azimuth, but takes into account all possible orbital 
inclinations. For eccentric orbits Pdetoct would be higher 
for any given a, as the planet would spend more time 
close to apastron where it is easier to detect. 

In Figure 4 we plot the detection probability for a 1 
Mj U p planet and a 2 Mj up planet using IWA = 4.5 AU 
for the 1 Mj U p case and IWA = 3 AU for the 2 Mj up case 
(see, lefthand plot in Figure 2). It is clear that we had 
a ~50% chance to detect (at greater than 5-<r) a young 
Jupiter analog at 5.2 AU, assuming one is present. The 
probability rises quickly for more massive or more distant 
planets. 

3.3. Comparison of theoretical models 

In the righthand plot of Figure 2, we have over-plotted 
in horizontal lines and shaded regions the magnitudes 
for a 1 Mj up object with an age of either 10 or 50 Myr 
based on theoretical predictions from the COND models 
(Baraffe et al. 2003) or more recent models published by 
Spiegel & Burrows (2012). The COND models predict 
that even for an age of 50 Myr we would have been able to 
detect objects with masses between 0.5 - 1 Mj up for sepa- 
rations >5 AU. The models by Spiegel & Burrows (2012), 
however, predict that planets are significantly fainter in 
L' compared to the COND models, regardless of their 
formation history (hot vs. cold start models). In this 
case, only for young ages (<20-30 Myr) are our detec- 
tion limits sufficient to detect a 1 Mj up object in the 
background limit. For an assumed age of 40 Myr only 
objects with masses >2 Mj up would have been detected. 

4. DISCUSSION 

Our detection limits are, to our knowledge, the deep- 
est yet achieved for a star targeted by direct imaging 
in terms of companion mass at a given physical separa- 
tion. Our results also illustrate an important challenge in 
interpreting (detections and) non-detections from high- 
contrast surveys: the systematic and strong dependence 
upon thermal evolutionary model predictions (see right 
panel in Figure 2). The Spiegel & Burrows (2012) mod- 
els consistently predict fainter L-band magnitudes for 
an 1 Mj up object compared to the Baraffe et al. (2003) 
COND models, most likely as a result of different as- 
sumptions regarding atmospheric opacities (Spiegel 2012, 
private communication). As mentioned in the introduc- 
tion, presently available atmospheric models fail to self- 
consistcntly explain the colors and luminosity of young, 
low surface gravity objects (e.g., Mohanty et al. 2007; 
Skemer et al. 2011; Barman et al. 2011; Skemer et al. 
2012). 

Nevertheless, our observations demonstrate that AP 
Col it is unlikely to host a gas giant at separations >5 
AU (Figure 4). For comparison, Delormc et al. (2011) 
have observed 14 nearby, young M-stars, each resulting 
in non-detections. We note, however, that the highest 
detection probability of that survey was achieved at sep- 
arations >30 AU. Monte Carlo simulations of current and 
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planned high-contrast surveys indeed indicate that low- 
mass stars are least likely to harbor a directly detectable 
planet (Crepp & Johnson 2011). These simulations are 
based on empirical results from radial velocity studies 
that find occurrence rates of <2% for giant planets with 
close-in orbits (Johnson et al. 2007; Cumming et al. 2008; 
Bonfils et al. 2011). 

However, microlensing results suggest that 17+g% of 
low-mass stars host a gas giant planet with mass be- 
tween 0.3 - 10 Mj up at separations between 0.5 and 10 
AU (Cassan et al. 2012). Thus, the planet frequency 
seems to rise between the orbital separations probed by 
Doppler observations and those by microlensing. Com- 
bining the currently available results from RV, microlens- 
ing, and direct imaging, Quanz et al. (2012) have shown 
that the population of gas giant planets with semi- 
major axes between ~0.03 - 30 AU can be described 
by cZ/pianot = M a a@dM da, with M being the planets' 
mass, a their semi-major axis, and a — —1.31 (cf. Cum- 
ming et al. 2008) and f3 < 0.5 — 0.6. To further constrain 
the value for /3 and the orbital separation range where it 
is valid, additional direct imaging observations effectively 
probing the regions between ~10 - 30 AU are required. 

5. SUMMARY & CONCLUSIONS 
Our findings can be summarized as the following: 

• We have used two different telescopes to obtain 
deep, high-contrast, L-band observations of the 
nearest (8.4 pc) pre-main sequence star, AP Col, 
to search for planetary mass companions. Neither 
dataset revealed the presence of a candidate com- 
panion. 

• Our derived detection limits are the most sensitive 
yet obtained by direct imaging in terms of planet 
mass for a given physical separation. For an as- 
sumed age of 40 Myr, we could have identified ob- 
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jects with masses > 2 Mj up at separations between 
3 - 4.5 AU, and > 0.5 - 1 Mj up for separations 
>4.5 AU using the predicted brightness of planets 
according to the COND models. 

• Accounting for the unknown inclination and posi- 
tion of a potential planet, we had a ~50% chance 
to directly image an 1 Mj up object at 5.2 AU, pro- 
vided a young Jupiter analogue exists. 

• Our detection limits depend on the predicted L- 
band magnitudes for young, planetary mass ob- 
jects. Using two different sets of theoretical model 
predictions, the limits change significantly, demon- 
strating that: (a) care need be taken for the inter- 
pretation of direct imaging results, and (b) empir- 
ical model constraints are urgently required. 

Dedicated imaging surveys for planetary mass com- 
panions around M-stars are ongoing (e.g., Dclorme et al. 
2011). They will eventually provide stringent constraints 
on the overall frequency and semimajor axis distribu- 
tion of gas giant planets at separations >10 AU, com- 
plementing the results of RV and microlensing studies. 
Our results support the notion that future high-contrast 
imaging programs would maximize their yield by pref- 
erentially selecting massive (A, F-type) stars as targets 
(Crepp & Johnson 2011). 
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Searching for young Jupiter Analogs around AP Col 
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Fig. 1. — Comparison between the 5— a contrast curves obtained with Keck/NIRC2 and VLT/NACO. 
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Fig. 2. — 5— a detection limit in L' magnitudes for our VLT/NACO data as a function of separation from AP Col (black crosses and 
curve). Left: Overplotted as horizontal dashed lines are predicted L' magnitudes for 40 Myr old objects with masses of 0.5 , 1 and 2 Mj up 
at the distance of AP Col based on the COND models (Baraffe et al. 2003). The red asterisk shows the location of a young Jupiter analog 
(i.e, 1 Mj up at 5.2 AU projected separation). Right: Overplotted are predictions from theoretical models for a 1 Mj up object with an age 
of 10 or 50 Myr. The blue, dashed lines are based on the COND models. The grey shaded areas are predictions from the hot and cold 
start models from Spiegel & Burrows (2012), where the spread in magnitude for a given model, i.e., the width of the grey shaded areas, 
shows the difference between atmospheres with and without hybrid clouds. All models assume solar metallicity. 



TABLE 1 

Summary of deep L' imaging observations of AP Col 



Parameter VLT/NACO Keck/NIRC2 

No. of detector reads X exp. time 64x0.2 s 300x0.106 s 

No. of data cubes / stacked images 134 cubes 67 images 

Parallactic angle start / end -50.11° / +59.70° -16.24° / +5.59° 

Airmass range 1.02 - 1.05 1.72 - 1.78 

Typical DIMM seeing 0.7"... 1.0" ~0.6" 

PSF FWHM a ~0.12" (0.10") ~0.09" (0.08") 

(£C) mC an b 19.1% 



a Measured FWHM and theoretical diffraction limited FWHM in parenthesis. 

b Average value of the coherent energy of the NACO/PSF in data cubes. Calculated by the Real Time Computer of the AO system. 
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Fig. 3. — Final VLT/NACO image after a fake "Jupiter" has been inserted in the raw images at two different positions with a projected 
separation of 5.2 AU from the central star. The brightness of the "planets" corresponds to a 1 Mj up object with an age of 40 Myr based on 
the COND models. Both "planets" are clearly detected with with a signal-to-noise of ~9.6 and ~12.1, respectively. The innermost regions 
that are dominated by subtraction residuals have been masked out. 
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Fig. 4. — Detection probability of a 1 Mj up object (black, solid line) and a 2 Mj up 
VLT/NACO observation assuming the detection limits from Figure 2 (left). The vertical, 
our Solar System at 5.2 AU. 
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